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ABSTRACT 

Peculiar velocities are one of the only probes of very large-scale mass density fluctua- 
tions in the nearby Universe. We present new "minimal variance" bulk flow measure- 
ments based upon the "First Amendment" compilation of 245 Type la supernovae 
(SNe) peculiar velocities and find a bulk flow of 249 ± 76 km s^^ in the direction 
I = 319° ± 18°, b = 7° ± 14°. The SNe bulk flow is consistent with the expectations 
of ACDM. However, it is also marginally c onsistent with the bulk flow of a l arger 
compilation of non-SNe pecuhar velocities (jWatkins. Feldman. fc Hudsonll2009f) . By 
comparing the SNe peculiar velocities to predictions of the IRAS Point Source Catalog 
Redshift survey (PSCz) galaxy density field, we find il'^^^as,\in = 0.40 ± 0.07, which 
is in agreement with ACDM. However, we also show that the PSCz density field fails 
to account for 150 ± 43 km s^^ of the SNe bulk motion. 



1 INTRODUCTION 



In the standard cosmological model, gravitational instability 
causes the growth of structure and peculiar velocities. In the 
regime where the perturbations are linear, there is a simple 
relati onship between density and peculiar velocity (| Peebles! 
1 19931 ): 



(1) 



where t he growth fa ctor / is equal to fim^^ in flat ACDM 
models (|Lindeill2005l l. 5 is the normalized mass density fluc- 
tuation field, 5 = [p — p)/p, and r are coordinates in units 
of km s~^. 

Given set of peculiar velocities, one can define a bulk 
flow as their average velocity; ideally the peculiar veloc- 
ity tracers are dense and numerous enough that the result- 
ing average is representative of the velocity of the volume. 
The bulk flow is then primarily due to structures on scales 
larger than the vol ume over which the bulk flow is measured 
(see Appendix A of lJuszkiewicz et al.ll 19901 for a derivation). 
Hence, bulk flows are probes of the large-scale power spec- 
trum of matter density fluctuations. 

The ACDM model, once normalized by WMAP7 
Ibarson et al.llioill ') observations of the Cosmic Microwave 
Background (CMB), fully specifles the r.m.s. fluctuations 
of 5 on all scales, an d hence the cosmic r.m.s. of bulk flows 
llWatkins et al .I2OO9I I. While most studies of bulk flows agree 



on the general direction of the flow, t here is some di s agree - 
ment as to the amplitude and scale. IWatkins et al.l (|2009l ) 
applied a "Minimal Variance" (MV) weighting scheme to a 
compilation of 4481 peculiar velocity measurements. Their 
results correspond to a sample with an effective Gaussian 
window of 50 Mpc and show a bulk flow of 407 ± 81 
km s"^ towards I = 287° ± 9°, & = 8° ± 6°, which is in con- 
flict with ACDM -I- WMAP7 at the 98 percent confldence 
level. The most controversial bulk flow re sult is the kinetic 
Sunyaev-Zeldovich flow dipole reported bv lKashlinskv et al.l 
(2010. ). who found a bulk flow on the order of 1000 km s~^ 
in the direction of I = 296° ± 28°, b = 39° ± 14° over a scale 
of at least 800 /i~^Mpc. If correct, this result would strongly 
conflict with ACDM + WMAP7. 

Another approach to understanding large-scale motions 
is to try to reconstruct the motion of the LG with re- 
spect to the CMB (627 ± 22 km s"^ towards I = 276° ± 3°, 
fe = 30° ± 2°: lKogut et al.ll 19931 ') by measuring the distribu- 
tion of galaxies and calculating the peculiar velocity of the 
Local Group (LG) using Eq. ([l]). Given the gravitational in- 
stability model of linear theory, the predicted velocity should 
converge to the measured CMB dipole for a sufficiently 
large survey volume. The application of Eq. ^ is difficult 
in practice because there are few r edshift surveys that are 
both a ll-sky and deep. For example. iRowan- Robinson et al.l 
(l200Cf l found that the predicted dipole from the IRAS Point 
Source Catalog Redshift survey (jSaunders et al.ll2000l . here- 
after PSCz) converged to 13.4 degrees of the CMB dipole by 
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30,000 km s"^ However. iBasilakos fc Plionid (|2006l ') reana- 
lyzed the same data set and found that significant power was 
required on large scales, which was missed by the original 
analysis. Other studies have been b ased on the Two Mi- 
cron All-Sky Sur vey Redshift Su r vey jSkrutskie et al.llioO^ . 
hereafter 2MRS ) : lErdogdu et~aD ('2006') and found probable 
convergence, but Lavaux et al.l {20 10) concluded that con- 
vergence was not obtained by 12,000 km s~^, and may not 
be until well beyond 20, 000 km s~^. In anot her study, using 
only the infrared fluxes, iBilicki et al.l l|201ll ') concluded that 
even at an effective distance ~ 300/i~^ Mpc (Ks < 13.5) the 
flux dipole had not converged. 

In this paper we use Type la SNe as our peculiar 
velocity tracers. SNe have also been us ed as peculiar ve- 
locity probes by a number of authors llRiess et ahl 1 19951 . 
1997i: iR adburn-Smith ct al.' '2004'; 'Pike & Hudson' boosl; 



Lucev et a l. 2005; HauEbocllo ct al. 2006; CoUn ct al. 2011; 
Dai et al.lboill ). 

An outline of this paper is as follows: In section 2, we 
introduce the data sets that were used. Section 3 presents the 
bulk flow of the SNe, using both simple weighting schemes 
as wel l as the "Minimal Variance" scheme of [Watkins ct al. 
l|2009l ). Section 4 compares individual SNe peculiar velocities 
to the predictions of the IRAS PSCz density field. We discuss 
the implications of our results in Section 5, and present our 
conclusions in Section 6. Throughout, we adopt Urn = 0.3, 
flA = 0.7, and quote distances in units of km s~^. 



we dub the 'First Amendment' (Al) compilation which we 
consider to be an extension to the 'Constitution' data selQ. 

Where the observed SNe in the data sets were known to 
be contained within a cluster of galaxies, the redshift of the 
cluster was used for the observed velocity distance rather 
than the redshift of the supernova itself Substituting clus- 
ter velocities for supernova velocities removes a significant 
source of thermal noise as objects in clusters can have a ve- 
locity rms of thousands of km s~^ . This process was applied 
to all three data sets. For galaxies not in clusters, the red- 
shift of the host galaxy was used if the host galaxy redshift 
was recorded in NED, which occurred in all but two cases. 
For the remaining two cases we used the redshifts of the 
SNe. Galactic longitudes and latitudes for the Carnegie set 
were also taken from NED. 

The Al data set has a characteristic or uncertainty- 
weighted depth of 58 h~^Mpc, where we define the charac- 
teristic depth to be: 



E 1/(72 



(2) 



where a is the total uncertainty in each SNe's peculiar ve- 
locity and r is the coordinates in units of km s~^. 

In Fig. [T] we present our results, our raw data, and the 
bulk flow directions that other surveys have found in an 
Aitoff projection. In Fig. [2] we present the Al data set in a 
Hubble Diagram divided into its three subsets. For all three 
data sets the intrinsic uncertainty of SNe's is the dominant 
source of error. Thus, for all our SNe the percent error is 
approximately 6 percent of the measured distance, with the 
scatter for the 'Old' and 'Hicken' subsets being larger. 



2 DATA AND CALIBRATION 

In this study, three primary data sets of nearby SNe (with 
distances less than 20,000 km s~^) are combined. 

We refer to the first of these data sets as the 'Old' sam- 
ple and it contains 106 SN the youngest bei ng from 2002, 



drawn from two sources: iJha et al.l l|2007l ) and lHicken eTal 
20091 Of the SNe in the 'Old' sample, 34 are fr omljha etld 



2007 ). The remaining 72 SN in 'Old' are from lHicken et al 



20091) . The second data set, which we refer to as 'Hicken' , 



contains the remaining 113 SNe from iHicken et all (|2009l ') 

after cutting objects at distances larger than 20,000 km s~^ 
and cutting two more objects (sn2007bz and sn2007ba) be- 
cause they deviated by more than 3a after the first round 
of fitting (as described below). The last set is the recently 
released data set fro m 'The Carnegie Supernova Project' 
llFolatelU et all I2OI0I . hereafter CSP), containing 28 SNe. 
Two of these objects were discarded due to our 20,000 
km s~^ distance cut, leaving 26 usable SNe. The CSP's re- 
ported uncertainties only refiected the derived distance mod- 
ulus residual spread. A second intrinsic uncertainty (agiv) in 
the magnitude of the SNe was added in quadrature by fit- 
ting a fiow model and reducing the reduced chi-squared fit 
to 1.00. The intrinsic uncertainty was found to be 0.107 mag 
(slightly smaller than the 0.12 mag found by the CSP due 
to cuts and the additional free parameters of bulk fiow). 
For further discussion of the light curve fitting, and conse- 
quences there of, for the 'Old' and 'Hicken' data sets see 
Appendix A. 

We combine these three sets to create a new sample that 



3 BULK FLOW 

In this section we discuss the bulk fiow, which is the simplest 
statistic that can be derived from a peculiar velocity survey. 



3.1 Methods 

We use two methods to measure the bulk flow. The first 
is a Maximum Likelihood (ML) method that minimizes the 
measurement uncertainties. The ML method is the tradi- 
tional method used and we consider it in order to compare 
new results with previous ML results. However, ML methods 
have the disadvantage of returning the bulk fiow of a spe- 
cific sparse sample of peculiar velocity tracers rather than 
the bulk fiow of a regular volume. Comparisons between 
ML results are complicated by the different spatial sam- 
pling. Instead, what is of greater interest is the bulk flow of 
a standardized volume. To estimate this, we calculate the 
"minimum var iance" (MV) bulk flow as flrst introduced by 
IWatkins et al] ((20091 



^ The 'Old' and 'Hicken' se ts combined r e sembl e very closely 
the 'Constitution' set from iHicken et al. in terms of 

which supernovae are included. The light curve fitter used here 
(MLCS2k2) differs from that of the 'Constitution' data set 
(SALT2). 



Cosmic flows in the nearby universe from Type la Supernovae 3 




Figure 1. An AitofT projection of our data with circles (asterisks) representing SNe witli peculiar velocities towards (away from) the LG. 
Larger symbols represent larger peculiar velocities in accordance with the scale shown top and bottom left. Also plotted in triangles are 
the direction motion of the LG with respect to the CMB, the Watkins et al. (2009) bulk flow direction, the kinetic Sunyaev-Zel'dovich 
bulk flow direction of lKashlinskv et alj | |2010| . labelled KAEEK), and our new results (labeled SNBF for the bulk flow results from Section 
3 and SNRF for the residual flow discussed in Section 4). 
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Figure 2. A Hubble diagram showing the three subsets that make 
up the Al data set: 'Old' (red filled circles) , 'Hicken' (Blue stars) 
and 'Carnegie'(green triangles). The error bars can be seen to be 
approximately constant in the log log diagram, or increasing pro- 
portionally to the distance, as to be expected with the dominant 
error in most cases being the intrinsic uncertainty in SNe. 



3.1.1 Maximum Likelihood 

In general, we fit a simple flow model (vpred) to tlie SNIa 
peculiar velocity data. In the case of the bulk flow V in the 
CMB frame, this flow model reduces to the radial component 
of the bulk flow vector for each SNe, i. e.: 

^^prod.i = V ■ f i (3) 

where rt is the unit vector pointing to each supernova. 

In the maximum likelihood method the weights are sim- 
ply determined by the total uncertainty on the peculiar ve- 
locity of each object. Uncertainties in the observed peculiar 
velocity can be approximated well by a Gaussian, in which 
case the Maximum-Likelihood solution can be approximated 
by the following x^- 

\cZof)s.i (^i ^^^prcd,^)] 



(4) 



where cZots is the observed redshift in km s^ , is the 
distance converted from the reported distance modulus in 
km s~^ , Uprcd.i is the model velocity we are trying to measure 
as predicted for SNe i, and ai is the total uncertainty on 
the peculiar velocity of object i in units of km s~^. This 
total uncertainty is the quadrature sum of the measurement 
error o"^,i, the intrinsic uncertainty on SNe magnitude asN 
(both converted from magnitudes to km s~^) and a "thermal 
noise" term (crth) in units of km s~^ due to uncertainties in 
the flow model, such that; 

f^i = '^m,i + (^SN + f^th (5) 

Note that since the ctIjv term is converted from an uncer- 
tainty on magnitude, it is proportional to the distance to 



4 S.J. Turnhull et al. 



the SN. 

(j|jY is often the dominant source of uncertainty 
since the thermal term is only important in nearby super- 
novae where the a'^^i and a|jv terms are small. The results 
are only weakly dependant upon the precise value chosen. 
Here, where the flow model is a simple bulk flow, we set the 
thermal noise to 250 km s~^, which is consistent with pre- 
vious work. The impact of this choice for the thermal noise 
is discussed in the results below. 

We let each sub-sample of the Al data set have a freely- 
varying independent Hubble term to identify degeneracies, 
to avoid underestimating final uncertainties, and to account 
for the fact that each subset may have slightly different cali- 
bration. None of the fits preferred a Hubble value that varied 
by more than 1 percent from the value drawn from the orig- 
inal sources. 

3.1.2 Minimum Variance Method 

While the ML method described above is the best estimator 
of the bulk flow of a sparse sample, it is restricted in that 
it can only really characterize a particular survey, that will 
have its own errors and a specific and somewhat ill-defined 
geometry. The ML method is also, in a sense, density sam- 
pled, with higher density regions being more likely to contain 
a SN than voids. Most importantly, because the weights in 
the ML method are determined by the uncertainty on posi- 
tion in km s~^, ML methods can be dominated by nearby 
SN that have smaller distance uncertainties. 

To better approximate a ^;o/ttme-weight ed bul k flow , 
we use the prescription described in Watkins et al.l (|2009l ) 
to estimate the volume flow. Each SN is weighted so as to 
minimize the variance between the bulk flow measured in 
the real sample and the bulk flow as it would be measured 
in a perfectly-sampled 3D Gaussian. We adopt a Gaussian 
with an "ideal" radius Ri = 50/i~^Mpc. Effectively, weights 
are assigned to each SN based on their proximity to other 
SNe in the data set, and on how they compare with an 
ideal uniform sampling. This weighting scheme is specifl- 
cally designed to maximize sensitivity to large scales. The 
MV weighting scheme has been tested using mock catalogues 
drawn from N-body simulations by Agarwal et al. 2011 (in 
preparation), who demonstrate that the recovered MV bulk 
flows are unbiased and have errors within the range expected 
from linear theory. 

3.2 Consistency of SNe Subsamples 

Before analysis of the combined SNe sample is undertaken, 
it is important to conflrm that the data subsamples agree 
with one another. We calculate a s tatistic for each pair o f 
subsamples, following the analysis of lWatkins et aL I (|2009D . 
which accounts for sparse sampling effects. The statistic 
we use is given by the equation, 

x' = ^(AV'0(a,,)-^(AV-). (6) 

where AT? is the bulk flow vector, and C is the covariance 
matrix taking into account the window functions of both 
su rveys and the power spectrum, (see in equations 21 - 23 
of lWatkins erail (|2009l ')'). The results are shown in Table [T] 
In summary, we flnd that all three subsamples are consistent 
with each other. 



Table 1. for 3 DoF for the surveys for fi™ = 0.258. If the 
value is greater than 7.8, the two surveys disagree at a greater 
than 95% confidence level. The probabilities resported are the 
test probability of agreement between the two. 



Ri = 


50/i 


-^Mpc 


Survey 




Probability 


Old vs Hicken 


0.173 


98.2 


Old vs Carnegie 


2.293 


51.4 


Hicken vs Carnegie 


1.369 


71.3 



3.3 Results 

In Table [2] we present a summary of the results from the 
bulk flow, subdivided by data set and by weighting scheme. 
The ML bulk flow for the Al sample was found to be 197 ± 
56 km s"^ in direction I = 295° ± 16°, 6 = 11° ± 14°. This 
is significantly different from zero at the 99.9% confidence 
level. 

As discussed above, the ML method gives most weight 
to SNe with the lowest errors in units of km s"'^, i.e. the 
nearest SNe. In order to reduce the impact of these nearby 
SNe, it is interesting to redetermine the bulk flow excluding 
nearby objects. The middle section of Table[2]shows the bulk 
flow using only SNe with 6000 km s~^ < d < 20000 km s"\ 
This subsample indicates a slightly higher amplitude flow, 
albeit with larger error bars: 330 ± 120 km s~^ towards 
I = 321°, b = 20°. 

Finally, the MV results shown in the third section of 
Table [2] should give the most robust estimates of the flow of 
a Gaussian volume of radius 50 /i~^Mpc. For the entire Al 
sample, the MV flow is 248 ± 87 km s^^ in the direction 
I = 319° ± 25°, 6 = 7° ± 13°. 

These values are lower than the LG's motion in the 
CMB frame, indicating that some of the LG's motion must 
come from structures within our survey volume (such as the 
Virgo and Hydra-Centaurus superclusters). 

To investigate the sensitivity of our results to the value 
of the thermal noise we adjusted it by ± 100 km s~^; When 
so tested, the final magnitude of the Al sample MV flow 
only changed by ± 31 km s~^. 



3.4 Bulk Flow: Cosmology and Comparisons 

It is in teresting to comp are our ML bulk flow result to 
that of IColin et al.l (|2011h , who apply a maximum likeli- 
hood bulk flow fit to th e Union2 catalogue of Type la SNe 
( Amarmll ah et al.ll2010^ . The Union2 catalogue contains 557 
SNe, of which 165 are within 30,000 km s"^ IColin et all 
(j201lh 's analysis yields a bulk flow velocity of 260 ±150 
km s~^ based on SNe within 18000 km s~^. Our Al sample 
yields a ML result of 196 ± 55 km s~^, which is consistent 
with theirs. It must be noted that the agreement between 
these results is not as significant as might at first be as- 
sumed because there is significant overlap between the data 
sets. However, Union2 uses SALT2 rather than MLCS2k2 
{Rv~1.7) to obtain SN distances from the light-curve data. 
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Table 2. Bu lk flow for all three SNIa data subsets and the combined First Amendment set. For comparison MV50 result from 
IWatkins et"al] (2009) arc also included. Note the uncertainties quoted for the ML method are the just propagated uncertainties from 
measurements, the uncertainties for the MV method also include an approximation of additional noise due to non-uniform sampling. 



Number 


Mag 


1° 


b° 




Vy 


Vz 


ML Thermal noise=250 km s^^ 








km s"-*- 


km s"-"- 


km s^-"- 


Old 106 


226±76 


307it21 


4il5 


136it76 


-180it81 


-35it59 


Hicken 113 


142±85 


283it41 


30ib37 


27-1-94 


-120it96 


71it78 


Carnegie 26 


260±140 


330±170 


76±38 


54±182 


-35±230 


250±150 


First Amendment 245 


196±55 


300±17 


15±14 


94±55 


-165±58 


50±44 


ML Thermal noise=250 km s'^, d > 6,000 km s 


-1 










Old 45 


450±190 


331±26 


6±21 


390±200 


-210±190 


44±160 


Hicken 76 


280±180 


313±33 


27±25 


170±170 


-180±190 


130±110 


Carnegie 15 


1132±850 


117±14 


16±20 


-490±540 


970±810 


310±300 


First Amendment 136 


330±120 


321±20 


16±15 


250±120 


-200±130 


94±84 


MV Weighting Rj = 50/i~^ Mpc Thermal noise= 


=250 km s- 


-1 








Old 113 


240±110 


318±26 


-4±21 


180±110 


-160±110 


-16±86 


Hicken 113 


250±110 


310±25 


5±20 


160±110 


-190±110 


20±85 


Carnegie 28 


250±150 


0±340 


81±43 


40±190 


0±240 


250±150 


First Amendment 254 


249±76 


319±18 


7±14 


186±75 


-162±77 


32±59 


Watkins et al. C2009') 4481 


407±81 


287±9 


8±6 


114±49 


-387±53 


57±37 



iDai et all l|201lf ) also analyzed the Union2 catalog, 
spliting it into two subsets. They defined a nearby set with 
132 SNe at z < 0.05 for which they found a bulk fiow of 
188t\ll km s"^ towards / = 290l|?°, b = 20tg° which also 
agrees well with our results. The remaining 425 high-z SNe 
show no significant bulk fiow. This is expected since the 
peculiar velocity errors are typically 6% of the distance to 
the source, and for this distant sample the errors per SNe 
measured in km s~^ are extremely large. 

Anot her interesting r ecent analysis of all peculiar veloc- 
ities is bv I Watkins et all (|2009l ). who studied peculiar veloc- 
ities mostly from Tully-Fisher, Fundamental Plane and SNe. 
They found that those subsamples had bulk flows consistent 
with each othei0. They combined the individual peculiar ve- 
locity samples into a "Composite" sample of 4481 peculiar 
velocity tracers, which was found to have a MV50 bulk fiow 
of 407 ± 81 km s"^ towards I = 287° ± 9°, 6 = 8° ± 6°. This 
result is inconsistent with ACDM at the 98% CL. However, 
their sample is not independent of ours. 103 of the 108 SNe 
which make up the "Old" subset of Al are common to both 
Al a nd Composit e , alth ough the latter takes SNe distances 
from iTonrv et al.) l|2003h . When all SNe data are removed 
from the 'Composite' data set, the two surveys become com- 
pletely independent, and can be compared using the same 
formal ism described i n Sect ion 3 of this paper and Section 
5.1 of IWatkins et all (120091 '). We find that the 'Composite 
excluding SNe' MV50 bulk fiow and the Al MV50 results 
are consistent with each other, although the agreement is 
marginal: of 6.4 for 3 directional degrees of freedom yields 
a 9 percent probability that the two results are consistent. 

Lastly, these results can be compared directly to the 



^ Except for the BCG sample of lLauer &: Postmanl l ll994l) . which 
was excluded from further analysis. 



expectations for a ACDM universe. A plot showing the ex- 
pectations of the one dimensional rms for perfectly sampled 
Gaussi an sphere can be found in the top thre e plots of figure 
5 from iFeldman. Watkins. fc HudsonI l|2010l ). The expecta- 
tion for the one dimensional rms for a perfect Gaussian is 80 
km s~^. When you take into account the sparse sampling of 
the real Al data set, this rises slightly to 91.2 km as- 
suming a as of 0.8 and flm of 0.258|j. If you then include 
the propagation of measurement uncertainties the total ex- 
pected rms for surveys equivalent to ours at different loca- 
tions ins space is 121 km s~^. This prediction leads to 
of 3.70 for 3 directional degrees of freedom yielding a 70 
percent probability that the Al data set is consistent with 
ACDM. 

4 PREDICTED GRAVITY FIELD 
4.1 Introduction 

The MV weighting scheme discussed above is designed to 
suppress the effects of small-scale flows that would other- 
wise "alias" power into the bulk-flow statistic. An alternate 
method for removing the effects of small-scale structure on 
flow measurements is to assume gravitational instability and 
linear perturbation theory Eq. ([1} and to predict the pecu- 
liar velocities using a model of the density fleld (derived from 
an all-sky galaxy redshift survey). The result is a model- 
dependent correction to measured peculiar velocities which 
can separate local effects from large-scale density waves from 
outside the survey volume. 

^ In the rest of the paper because SN distances are insensitive to 
the value of Qm we used Qrn, = 0.3. The predicted one dimensional 
rms drops slightly to 88.4 if this slightly higher value of Qm is used 
in the prediction 
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Suppose we have an all-sky redshift survey that extends 
to a distance -Rmax- We will model the peculiar velocity of a 
given SN located at position r by setting iiprcd.i of Eq. Q 
to a function with two terms: 



SGZ = km plane 



/3 



Vp.ed(r) = — 



d\'5,{v') 



+ U 



(7) 



where (3 = f /b, fe is the linear bias between galaxy density 
and mass density, U is the residual bulk flow of the vol- 
ume being driven by mass structure beyond -Rmax. In prin- 
ciple, the re sidual velocities have tidal (shear) and higher 
order terms. iFeldman et"all (|2010l ) measured the the tidal 
and higher order te rms for the "Composite" sample of 
IWatkins etlll l|2009l '). but found them to be small. We ne- 
glect these terms here and model the residual as a simple 
bulk flow U. 

The first term of Uprod is the predicted peculiar veloc- 
ity induced by structure within the redshift survey volume 
(r < J?max). The model is scaled by P to match the ob- 
served peculiar velocities of the SNe tracers. The peculiar 
velocity data therefore yields information about Q,m and 6. 
The residual bulk flow U is the additional velocity of the 
entire redshift survey volume in the CMB reference frame, 
and is presumably due to sources beyond -Rmax- In an ideal 
survey, U would be completely independent of any structure 
within -Rmax. This de-coupling of U from /3 means that U 
can be used to test consistency with ACDM + WMAP7 on 
large scales and /? can do so on smaller scales. 



4.2 Data and Method 

The PSCz is both all-sky and d eeper than, for example, 
the 2MRS ijSkrutskie et all bOOd). Here we us e the PSCz 
density field reconstructed bv ISranchini et "all () 19991 ). For 
our study we applied the same 20,000 km s~^ limit to the 
PSCz as we applied to the SNe. The PSCz density field in 
the Supergalactic Plane is shown in Fig. [S] 

We fit the SNe data using the same method as in Section 
13.1.11 but now with a new model as given by Eq. ((?]). Since 
the integral is specified by the PSCz density field, the free 
parameters are /3 and the three components of U. Since the 
PSCz plus bulk fiow is a better fiow model than a simple bulk 
fiow, we reduce the thermal compon ent to 150 km s~^, w hich 
is consistent with previous studies (|Hudson et al.|[2004l ). 



4.3 Results 

The results of the fits to each subset are given in Table[3l We 
find that the results from independent subsets are consistent 
with each other. For the Al sample, the magnitude of the 
residual bulk fiow was found to be 150 ± 43 km s~^ in 
direction I = 345° ±20°, 6 = 8° ± 13°. This is significantly 
different from at the 99.6% CL. 

The value of /3 was found to be 0.53 ± 0.08, and is shown 
in Fig. U The fit is sensitive to a single outlier, snl992bh, 
for which the PSCz prediction is rather high (1719 km/s). 
Excluding this SN, we find P = 0.57 ± 0.08. 

Again we investigated the sensitivity of these results on 
the thermal noise term by changing it by ± 100 km s~^. 
Again the magnitude off the fiow only changed by ± 20 
km s~^, and the /3 changed by ± 0.03. 




SGX (km s 



Figure 3. The Supergalactic Plane. The PSCz galaxy density 
field is shown by the contours, predicted peculiar velocities as 
small black arrows, and measured supernova positions as "tad- 
poles" with dots showing measured positions and tails showing 
the magnitude of the measured radial peculiar velocity. The thick 
black contour corresponds to a <5 = (or contours where the den- 
sity is the mean universal density). The rcd(filled circles) SNe 
have peculiar velocities away from the LG and the blue(open cir- 
cles) SNe have peculiar velocities towards the LG. 



Representative 




Vpred (km/s) 

Figure 4. The observed peculiar velocity minus the measured 
bulk flow as a function of the linear-theory-predicted peculiar ve- 
locity for each SN, assuming /3 = 1. The circular symbol diameter 
scales with the inverse of the uncertainty (hence symbol area is 
proportional to weight). Representative error bars are shown in 
the top left. The slope is the fitted /3 = 0.53. 
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Table 3. Results from all three data subsets and the Al full set with a linear perturbation theory model fit with a known matter 
distribution to fit /3 and the residual flow U. Fit with /3 as a free parameter and with a thermal noise of 150 km s-l. Note the 
uncertainties quoted for the this method arc the just propagated uncertainties from the measurements. 





# of SN 


Mag 


1° 


b° 




Uy 


Uz 


/3 




Dof 












km s"-*- 


km s"-*- 


km s"-*- 








Old 


106 


190±59 


349±22 


0±14 


187±60 


-36±73 


0±46 


.45±.ll 


139 


101 


Hicken 


113 


86±77 


347±54 


9±41 


84±87 


-19±83 


13±66 


.62±.13 


102 


108 


Carnegie 


26 


290±150 


347±41 


31±26 


240±170 


-50±190 


151±130 


.82±.33 


21 


21 


First Amendment 


245 


150±43 


345±20 


8±13 


144±44 


-38±51 


20±35 


.53±.08 


270 


238 



4.4 Gravity field: cosmology and comparisons 

As noted above, the residual bulk flow U is significantly 
different from zero at tiie 99.6% confidence level. This means 
that there are structures not found in the PSCz catalogue 
that contribute significantly to t he total peculiar velo city of 
the LG. As discussed in detail in [Hudson et al.l (|2004l ) these 
structures could be missing from the PSCz because they are 
(i) outside the survey volume (ii) in the Zone Of Avoidance, 
or (iii) present but underrepresented. The latter scenario 
may arise because the IRAS (far-infrared) selection on which 
PSCz is based is sensitive to d usty spirals, but less so t o 
the mostly dust-free early types. iLavaux fc HudsonI l|201ir ). 
using the 6dF survey of 2MASS-selected galaxies, showed 
that the Shapley and Horologium-Reticulum superclusters 
generate significantly more peculiar velocity than predicted 
by the PSCz, even allowing for a different /3 for 2MASS 
galaxies. 

We found that, for IRAS-selected galaxies, /3 = 0.53 ± 
0.08. This value of (3i is in good agree ment with the 
IRAS average of 0.50 ± 0.02 reported by IPike fc HudsonI 
ll2005ll.and with the SNIa-b ased resuh /3i = 0.55 ± 0.06 of 
iRadburn-Smith et al.l (|2004l l. Comparing fitted values of /? 
between redshift surveys of different galaxy types is com- 
plicated by the fact that the bias factor b need not be the 
same because different galaxy types may trace the under- 
lining mass density differently. This problem can be allevi- 
ated by noting that in linear theory the r.m.s. fiuctuation 
of the survey galaxies, say in an 8 Mpc top hat sphere 
(o"8,gai) is proportional to the true matter r.m.s. fiuctua- 
tions in a volume of the same size (i.e. (Tg.gai = berg, mass). 
Thus with our measured / jj and the known asj from the 
IRAS PSCz of 0.80 ± 0.05 (|Hamihon fc Tegmark|[20o3 ). we 
can calculate the degenerate parameter pair fas (where we 
dropped the subscript 'mass'). Our value of (3 corresponds 
to fas = 0.424 ± 0.069. 

We can the n compare our fas to other studies. 
iDavis et all(|201ll ) compared the 2MRS density field and the 
SFI-I— I- peculiar velocity data, and derived fas = 0.31±0.05. 
This is lower than our result, but not significantly so {1.5a). 

The fas parameter can also be derived from WMAP7 
results. Recall that WMAP is observing fiuctuations at an 
early epoch, when the perturbations were still well in the 
linear regime. To compare to WMAP7, we can convert our 
non-linear as into the equivalent lin ear value using the pre- 
scription of lJuszkiewicz et al.l l|2010l ). If we assume an of 
0.272, our as,un becomes 0.814 compared to its non-linear 
value of 0.867. Using this value of (Tg.un, /fg.iin drops to 



0.40 ± 0.07, which is in excellen t agreement with the results 
of WMAP7 l|Larson et al.ll201ll ): fas,iin = 0.39 ± 0.04. 



5 DISCUSSION 

Attempts to determine the sources of the LG's motion 
amount to determining the factors in Eq. ([l]). While early 
studies focussed on simple toy infall models, more recent 
studies have concentrated on models of the density field with 
the two free parameters /3 and U. For a single object, such as 
the LG itself, there is a trade-off between these parameters. 
Lower values of /3 lead to larger values of U, which are re- 
quired in order to match the same v on the left-hand side of 
equation [T] This degeneracy can be broken with more than 
one measurement. We have shown that the PCSz does not 
account for all of the motion of the LG, although it is plau- 
sible that some of the missing signal comes from within the 
survey volume in the form of extra infall into the highest- 
density superclusters. 

An alternative explanation for the bulk fiow has been 
proposed, namely that the CMB temperature dipole, or 
part thereof, is intrinsic and does not repres ent the pecu- 
liar velocity of the LG (A "tihed" Uni verse: iTurnerl Il99ll : 
iKashlinskv et aL I l2008l : iMa et all [20T1I '). This would lead 
to an illusory "bulk fiow" which would extend well be- 
yond the local volume, indeed to the horizon. The apparent 
1005 ±267 km s~ ^ bulk fiow of 2 < 0.25 clusters claimed by 
iKashhnskv et all (|2010l ). which is well outside the expecta- 
tions of ACDM bulk fiows, might be explained by such an 
effect. In such a scenario, there is an additional "bulk fiow" 
Utiit which never vanishes no matter how deep a redshift 
survey -Rmax is used in Eq. 0. Our measured U thus pro- 
vides an upper l imit on the Utiit ■ The amplitude of the bulk 
fiow found by i Kashlinskv et al. 1 Ikoio) is inconsistent with 
our measurem ent of U = 150 ± 43 km s~^ . However, am- 
plitude of the IKashlinskv et all (|201G| ) bulk flow is system- 
atically uncertain. If we compare only th e direction of the 
Al fit Z = 345° ± 20°, fe = 8° ± 13°and the lKashhnskv et al.l 
(2010) direction I = 296° ±29°, b = 39° ± 15°, the resuhs 
are marginal: they disagree at approximately the 90% CL. 
Thus our res ults do not suppo r t the high amplitude bulk 
fiow found bv IKashlinskv et all (|201G| '). 



6 CONCLUSION 

We have analyzed the peculiar velocities of a 245 SNe 
dataset dubbed the "First Amendment". Overall, we have 
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found that this new compilation is in marginal agreement 
with previous bulk flow results and is not in significant con- 
flict with ACDM + WMAP7 predictions. The First Amend- 
ment compilation yields a bulk flow of 248 ± 87 km s~^ in 
the direction / = 319° ± 25°, fe = 7° ± 13°. 

We have compared the peculiar velocities to the pre- 
dictions from the IRAS PSCz and have found fJm^^crg.iin of 
0.40 ±0.07, which is in excellent agreement with the ACDM 
-I- WMAP7 predictions and other previous measurements. 

A residual flow of 150 ± 43 km s'^ I = 345° ± 20°, 
6 = 8° ± 13° was found for the IRAS Point Source Catalog 
as normalized with the First Amendment SNe. This may 
suggest that the IRAS PSCz undersamples massive dense 
superclusters such as the Shapley Concentration. Neverthe- 
less, the small amplitude of the residual flow is in conflict 
with "tilted Universe" scenario s such as might be favoured 
by the kSZ analysis of Kashlins kv et al.l (|2010l ). 

As its name suggests, the First Amendment compi- 
lation is readily extendible as new SNe are found and 
their distances are published. Ongoing surveys such as 
CfA4 (95 SNe, Hick e n, pr ivate communication), LOSS 
( Ganeshalingam et al.l 1201 ij). Palomar Transit Factory 
( Law et aLll2009l '). and CSP l|Stritzinger et al.1l201ll , 50 more 
distances expected soon), and upcoming sur veys such as 
SkyM apper (100 SNe per year with z < 0.085, iKeller et all 
[2003), Pan-Starrs, and LSST will eventually provide suf- 
ficient SNe to reduce the 20,000 km s~^ bulk and residual 
fiow uncertainties to the systematic limits. Future results on 
/o"8,iin are expected based on predicted peculiar velocities 
from the 2M-I— I- redshift compilation (iLavaux fc HudsonI 
[ioij). Additionally, although individually less precise. Fun- 
damental Plane distances and peculiar velocities can con- 
tribute significant precision to bulk flow surveys by sheer 
numbers . We wish to r e-analyze the full 'Composite' data 
set from iFeldman fc~W atkins (200l) after replacing the 103 
SNe currently contained in that data set with the 245 SNe of 
Al, as well as to a dd the Fundament al Plane pe culiar veloc- 
ities from NFPS JSmith et al.ll2006l ) and 6dF (|jones et all 
l2004f ) when they become available. For now the results are 
data-limited, but the future promises many fruitful results 
from many promising surveys, and we await them eagerly. 
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APPENDIX A: LIGHT CURVE PARAMETER 
COMPARISONS 

The Al data set is composed of three different SNe cata- 
logues; this complicates the description of the light curve 
fitting procedures used because the catalogues u sed differ- 
ent rn ethods. In the 'Old' sample are 34 SNe from ljha et alJ 
|2003); most of which are fit using the MLCS2k2 light curve 
fitter with a reddening law parameter Rv of 3.1 (for SN 
with high extinction, Rv was a free fit parameter with a 
tight prior of 3.1). The remaining 72 SN in 'Old' are from 
iHicken et aD (j2009t ) and are also fit using MLCS2k2, but 
with a reddening law parameter Ry of 1.7. The second data 
set, whic h we refer to as 'Hick en', contains the remaining 113 
SNe from lHicken et all (|2009l ) and are all fit with MLCS2k2 
with a reddening law parameter Rv of 1.7. The 'Carnegie' 
set containing 28 SNe were fit with a Rv as a free variable. 
The light curve fitter used for the 'Carnegie' set is de scribed 
in detail in the ori ginal paper ISaunders etlo] (|2000l '). 

For SNe fit bv faicken et all l|2009l 'ir four distances were 



reported for each SN. We use the distances reported us- 
ing the MLCS2k2 fitting procedure rather than either of the 
SALT procedures for multiple reason. To start the MLCS2k2 
process determines host reddening on a case-by-case basis. 
Furthermore, of the two published MLCS2k2 methods we 
use the results with a reddening law p arameter Rv of 1.7 



instead of 3.1 since iHicken et ah (2009) show that the Hub 



ble residuals for high-extinction SNla's using Ry—S.l are 
systematically negative, (suggesting that the extinction is 
overestimated). We study the effect the choice of Rv pa- 
rameter has upon bulk flow measurements to explore sys- 
tematics. IHicken et al] l|2009l ) provide distances to 162 SNe 
using both Rv = 1.7 and Rv = 3.1. We fit both of these 
data sets for bulk fiows using the ML method to investigate 
the systematics. The results of this comparison can be seen 
in table lAll Although the results for the two light curve 
fitters agree to less than one a in each of the 3 degrees of 
freedom, the data sets are fit to the same light curves, so 
they are not independent. This result highlights how large 
the systematic errors are for bulk fiow surveys, in part re- 
flected by the large asN, which in most cases dominates the 
uncertainty budget for peculiar velocity surveys. 
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Table Al. Results for 162 SNe from lHicken et al.l l|2009l l fit with the MLCS2k2 Hght curve fitter either with i?v = 1-7 or i?v = 3.1. 



Number Mag 


1° 


b° 


Vx 


Vy 




ML Thermal noise=250 km s^^ 






km s^-*- 


km s^-"- 


km s^^ 


Rv = 1.7 162 220±70 
Rv = 3.1 162 175±70 


298±18 
310±25 


9±14 
14±18 


103±68 
108±70 


-191±73 
-131±75 


35±52 
43±53 



